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PREFACE 


Meteorological  studies  suggest  that  technologically  feasible  oper¬ 
ations  Bight  trigger  substantial  changes  In  the  cllaate  over  broad  re¬ 
gions  of  the  globe.  Depending  on  their  character,  location,  and  scale, 
these  changes  Bight  be  both  deleterious  and  Irreversible.  If  a  foreign 
power  were  to  bring  about  such  perturbations  either  overtly  or  covertly, 
either  aallclously  or  heedlessly,  the  results  might  be  seriously  detrl- 
mental  to  the  security  and  welfare  of  this  country.  So  that  the  United 
States  aay  react  rationally  and  effectively  to  any  such  sctlons.  It  Is 
essential  that  we  have  the  capability  to:  (1)  evaluate  all  consequen¬ 
ces  of  a  variety  of  possible  actions  that  Bight  modify  the  climate, 

(2)  detect  trends  In  the  global  circulation  that  presage  changes  In 
the  climate,  either  natural  or  artificial,  and  (3)  determine.  If  pos¬ 
sible,  aeans  to  counter  potentially  deleterious  cllaatlc  changes.  Our 
possession  of  this  capability  would  sake  incautious  experimentation 
unnecessary,  and  would  tend  to  deter  Bslicious  manipulation.  To  this 
end,  the  Advanced  Research  Projects  Agency  initiated  a  study  of  the 
dynamics  of  climate  to  evaluate  the  effect  on  cllaate  of  envlrons>sntal 
perturbations.  The  present  Memorandum  is  a  technical  contribution  to 
this  larger  study. 

Rand's  position  on  climate  and  weather  aedlf ication  studies  was 
asserted  in  its  publications  RM-3205-NSF  and  RM-5833-NSF.  The  approach 
to  understanding  the  consequences  of  cllaate  change  must  consist  of 
many  converging  paths;  and  for  reasons  of  safety  and  economy  this  ap¬ 
proach  must  take  advantage  of  the  growing  versatility  of  computers. 

This  Memorandum  is  the  third  in  a  series  reporting  investigations 
of  the  wind-driven  oceanic  circulation,  an  understanding  of  which  is 
essential  to  broadening  our  knowledge  of  the  ocean /atmosphere  inter¬ 
actions  that  so  powerfully  influence  weather  and  climate.  The  work  is 
based  on  models  for  solution  on  high-speed  computers.  Rand's  prior 
publications  on  this  particular  topic  are  RM-6110-RC  and  RM-6210-ARPA. 
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SUMMARY 


The  time-dependent  primitive  equations  for  e  shallow  homogeneous 
oceen  with  e  free  surface  ere  solved  for  e  bounded  basin  on  the  sphere, 
driven  by  e  steedy  zonal  wind  stress  and  subject  to  leteral  viscous 
dissipation.  These  are  the  vertically  lntegreted  equations  for  a  f ree- 
surfece  model,  and  ere  integrated  to  60  days  from  an  Initial  state  of 
rest  by  en  explicit  centered-difference  method  with  zero-slip  leterel 
boundery  conditions.  In  a  series  of  comperatlve  numer ice 1  solutions 
It  is  shown  that  at  least  a  2 -deg  resolution  is  needed  to  resolve  the 
western  boundary  currents  edequately  end  to  evold  undue  distortion  of 
the  transient  Rossby  weves.  The  8-plane  formuletlon  Is  shown  to  be 
an  edequate  approximation  for  the  mean  circulation  In  the  lower  and 
middle  latitudes,  but  noticeably  intensifies  the  transports  poleward 
of  about  50  deg  and  both  slows  end  distorts  the  transients  In  the  cen¬ 
tral  basin.  The  Influence  of  the  (southern)  zonal  boundary  on  the 
transport  solutions  Is  confined  to  the  southernmost  gyre,  except  In 
the  region  of  the  western  boundery  currents  where  Its  Influence  spreads 
to  the  northern  edge  of  the  basin  by  30  days.  The  total  boundary  cur¬ 
rent  transport  Is  shown  to  be  approximately  proportional  to  the  zonsl 
width  of  the  basin  and  Independent  of  the  basin's  (uniform)  depth, 
while  the  elevation  of  the  free  water  surface  Is  Inversely  proportional 
to  the  basin  depth,  In  accordance  with  llnaer  theory.  The  Introduction 
of  bottom  friction  has  a  marked  damping  effect  on  the  trenslent  Rossby 
waves,  end  also  reduces  the  maximum  boundary-current  transport.  The 
solutions  throughout  are  approximately  geos trophic  and  are  only  slight¬ 
ly  affected  by  the  nonlinear  Inertial  terms  for  simple  assumed  vertical 
profiles  of  the  horizontal  current. 

The  root-mean-square  (rms)  transport  verlabllity  during  the  period 
30  to  60  days  is  concentrated  In  the  southwest  portion  of  the  basin 
through  the  reflection  of  the  trenslent  Rossby  weves  from  the  western 
shore  and  has  a  maximum  corresponding  to  en  rms  current  variability  of 
about  3  cm  sec  The  transport  variabilities  are  about  10  percent  of 
the  mean  zonal  trensport  end  more  than  100  percent  of  the  mean  meridi¬ 
onal  transport  over  a  considerable  region  of  the  western  basin  (outside 
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the  western  boundary  current  regime) .  Some  99  percent  of  the  total  ki¬ 
netic  energy  Is  associated  with  the  zonal  mean  and  standing  zonal  waves, 
which  are  also  responsible  for  the  bulk  of  the  meridional  transport  of 
zonal  angular  momentum.  Although  the  transient  Rossby  waves  systemat¬ 
ically  produce  a  momentum  flux  convergence  at  the  latitude  of  the  maxi¬ 
mum  eastward  current,  much  In  the  manner  of  their  atmospheric  counter¬ 
parts,  this  Is  only  a  relatively  small  contribution  to  the  zonal  oceanic 
momentum  balance;  the  bulk  of  the  mean  zonal  stress  Is  here  balanced  by 
a  nearly  stationary  net  pressure  torque  exerted  against  the  meridional 
boundaries  by  the  wind-raised  water.  In  an  ocean  without  such  boun¬ 
daries  the  role  of  the  transient  circulations  may  be  somewhat  more  Im¬ 
portant. 
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I.  INTRODUCTION 


In  previous  studies  (Gates,  1968,  1969)  the  transverse  Rossby  wave 
was  shown  to  be  the  principal  node  of  response  of  a  homogeneous  ocean 
to  an  imposed  large-scale  wind  stress,  and  the  development  of  the  steady 
circulations  was  examined  with  emphasis  on  the  reflection  of  thesa 
transients  from  the  basin's  walls.  These  integrations  were  mada  for 
an  ocean  of  uniform  depth  on  the  B-plane,  with  a  modest  degree  of  non¬ 
linearity  and  a  prescrlbad  lateral  eddy  viscosity.  The  role  of  tha 
nonlinear  inertial  terms  and  the  effects  of  changes  in  the  viscous  dis¬ 
sipation  in  such  a  model  hava  baen  examined  by  Bryan  (1963),  whila  tha 
corresponding  effects  in  a  6 -plane  model  with  bottom  friction  hava  baen 
examined  by  Veronis  (1966c) .  It  is  the  purpose  of  the  present  paper 
to  extend  these  investigations  to  the  case  of  a  bounded  ocean  basin  on 
the  spherical  earth,  with  particular  attention  being  given  to  the  role 
of  the  transient  planetary  circulations.  In  addition  to  examining  the 
frictional  and  nonlinear  affects  in  the  spherical  case,  a  series  of 
comparative  numerical  Integrations  will  be  presented  to  show  the  effects 
of  resolution  (grid  mesh  size)  changes,  and  of  changes  of  the  basin's 
dimensions.  The  present  studies  of  a  homogeneous  ocean  may  be  vlewad 
as  part  of  a  continuing  series  of  experiments  which  will  culminate  in 
the  (barotroplc)  simulation  of  the  circulation  of  an  ocean  basin  with 
realistic  lateral  and  bottom  geometry.  Although  there  may  be  important 
barocllnlc  effects  in  the  large-scale  ocean  circulation,  the  barotroplc 
model  may  be  able  to  account  for  most  of  the  variability  of  large-scale 
currents  (Phillips,  1966a),  at  least  in  the  open  sea  of  middle  lati¬ 
tudes  . 
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II.  THE  DYNAMICAL  MODEL 


The  dynamical  bases  of  the  barotroplc  model  are  the  horizontal 
equations  of  notion  and  the  equation  of  continuity  for  an  incompres¬ 
sible  homogeneous  ocean.  In  the  spherical  coordinates  X  (longitude), 
$  (latitude),  and  z  (upward)  these  may  be  written 


3u  u  3u  ,  v  3u  ,  jhi 

3t  a  cos  4  3X  a  3$  W  3z 


uv 


tan 


4  - 


«  P. 


— - -  +  fv  +  — 

cos  $  3X  p 


3tX  2 

t— ^  +  A  V  u 
3z  s 


(1) 


2*  + 
3t 


u  3v  .  v  3v  .  3v  .  u 

3X  W  3l  +  ^tan  ♦ 


— ~  If  “  fu  +  —  T-^  +  A  V2v 
a  pQ  34  PQ  3z  s 


where  the  velocity  components  are  given  by 


(2) 


a  cos 

A  — 

9  dt 

(3) 

«£i 

*  dt 

(4) 

ds 

dt 

(5) 

with  a  the  earth's  radius,  pq  the  ocean's  (uniform)  density,  p  the 
pressure,  f  -  2D  sin  4  the  Coriolis  paraneter  (with  ft  the  earth's  angu¬ 
lar  speed  of  rotation) ,  and  the  eastward  and  northward  components 

of  the  tangential  stress,  A  the  coefficient  of  lateral  eddy  viscosity, 

2 

and  7^  the  spherlcml  Laplaclan.  Here  we  have  neglected  the  terms 

uw/a  *  and  2ftw  cos  4  in  (1)  and  the  term  vw/a~*  in  (2)  in  order  to 
maintain  energetic  consistency  when  the  hydrostatic  equation 
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“V 


(6) 


la  uaed  aa  tha  reanant  of  tha  vortical  aquation  of  action.  Thaaa  ap- 
proxlaatlona ,  togathar  with  tha  uaa  of  (a)  rathar  than  (a  ♦  a)  aa  a  co- 
afflclant  In  tha  taraa  of  (1)  and  (2),  hava  been  ahovn  to  ba  cona latent 
with  tha  praaarvatlon  of  angular  nouantun  (Phllllpe,  1966b).  Tha  aqua¬ 
tion  of  continuity  accompanying  (1)  and  (2)  may  ba  written 


1  »u  .  1  )  ,  .  v  .  >w 

HJT7  n  *  rST*  jj  <»«•♦)♦  •  « 


<»> 


whldi  axpraeaee  tha  nondlvorganca  of  tha  thraa-dlmanalonal  velocity. 
Hara  tha  extremely  nail  tarn  2wa  1  haa  baan  neglected.  Uelag  th la 
aquation,  tha  advectlve  taraa  (tha  aacond,  third,  and  fourth  on  tha 
laft-hand  aldaa)  of  (1)  and  (2)  nay  ba  rewritten  la  tha  "flux  forma" 
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Wa  now  latagrata  tha  ayatan  (1),  (2),  and  (7)  over  the  total 
vertical  depth  of  tha  hydroatatlc  ocean,  i.a.,  from  tha  aa tuned  rigid 
botton  at  a  •  -h(X,#),  at  which  tha  boundary  condition  la 


w-h 


9  h 

a 


(10) 


to  tha  free  upper  aurfaco  at  a  •  ((X,*,t),  at  which  tha  boundary  con¬ 
dition  la 


w  •  V 

C  Jt  \ 


V 


(11) 


Here  v  .  dtOktti  v  it  i  ■  -h,  V  .  li  the  horizontal  velocity  at  i  •  *h 
•n  -fl 

(vlth  cor  ii^oodl'4  lntarpratatlona  at  tha  surface  s  •  ()  ,  and  la 
the  spherical  gradient  operator.  Using  these  conditions,  the  hydro- 
static  equation,  (4),  and  the  aastaptlun  of  a  uniform  surface  (atmos¬ 
pheric)  pressure  et  s  •  (,  and  paying  particular  attention  to  the  vari¬ 
ability  of  the  limits  of  the  vertical  Integration  operator,  defined  es 


<(  » 


((>,♦. t) 


-**<».♦) 


(  )  d* 


(12) 


the  reeults  of  the  Integration  of  the  **flu*  forms**  of  (1)  and  (2)  ere 
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Here  t  and  r  denote  the  svrfece-erlnd  stress  and  the  tangential  stress 
on  the  oceen  bottom,  respectively.  A  slmller  Integration  of  the  con¬ 
tinuity  equation,  (7),  gives 
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1 

a  cos  4 


(15) 


Ths  furthsr  rsduccloo  of  chs  system  (13)  Co  (15)  hinges  upon  Che 
Introduction  of  chs  assumption  chsc  chs  lncsgrsls  of  Chs  quadratic  pro¬ 
duces  u  ,  uv  snd  v  occurring  In  chs  insreial  Carats  on  chs  loft-hand 
sldas  of  (13)  and  (14)  may  ba  expressed  in  carms  of  cha  products  of  chs 
lncsgrsls  (u)  and  (»);  this  assumption  amounts  eo  a  specification  of 
cha  vertical  structure  of  cha  horltontal  currant.  Specifically,  we 
shall  assuma 


<u2)  -  f  <u)2 
<uv)  •  ^  <u)(v) 

<v2>  -  £  <v>2  (16) 


where  S  Is  a  "shape  factor,"  assumed  to  ba  a  dimens  ionises  constant 
characteristic  of  the  ocaan,  and  h  is  tha  (undisturbed)  local  ocaan 
depth.  Tha  relations  (16)  ara  equivalent  to  assuming  that  tha  veloc¬ 
ity  components  vary  with  t  in  a  similar  fashion  in  every  vertical,  but 
at  a  rata  inversely  proportional  to  tha  local  watar  depth,  l.e., 

u  •  £  <u)  ,  v  -  £  <v)  (17) 


where  C  Is  a  dimensionless  function  of  s  with  (c)  ■  h.  Ths  shape  fac¬ 
tor  S  may  tharaforc  ba  written 
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If  a  currant  la  assumed  to  dacraaaa  linearly  toward  aero  at  the  ocean 
bottom,  for  example,  we  have  C  -  2(1  +  x/h)  and  S  -  4/3,  which  may  be 
compared  with  the  caae  of  a  current  unlforn  over  depth  for  which 
C  -  S  •  1.  For  a  parabolic  dacraaaa  of  apeed  with  depth,  with  a  poa- 
elble  reveraal  of  direction  in  the  deeper  water  (C  <  0) ,  valuea  of  S 
from  around  2  to  5  may  be  Inferred. 

Ulth  the  approximation  (16)  and  (18)  we  may  write  the  Integrated 
equation  of  motion,  (13)  and  (14),  in  the  forma 
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(20) 


Here  we  have  dropped  the  bracketa  and  written  u  and  v  for  (u)  and  (v) 

In  the  lntereat  of  simplicity,  and  we  have  Introduced  the  approximation 
C  «  h  in  the  coefficient  of  the  eecond  (pressure-force)  terms  on  ths 
right-hand  sides.  He  have  also  aeaumed  that  ths  last  (lateral  eddy 
viscosity)  terms  may  be  written  with  sufficient  accuracy  by  frssly  com¬ 
muting  the  3/3X,  3/34,  and  (  )  operators.  Togsthsr  with  the  continuity 
equation,  (15) ,  now  written  as 
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1L  1  /  3u  3(v  cos  i)  \ 

at  "  ‘  •  co«  ♦ v  ax  a#  / 


(21) 


thas*  aquation*  now  conatltut*  a  dynamical  ay a tan  in  the  volume  trana- 
port  coaponenta  u  and  v  (vertically-integrated  apeeda)  and  in  the  dia- 
placemant  of  the  free  ocean  aurface  (.  This  ayatea  la  the  dynamical 
baala  of  the  preaent  barotroplc  model,  although  it  haa  been  aimpllfled 
In  the  Integration*  preaented  below  by  the  aaaumptlon  of  a  uniform 
ocaan  dapth,  by  the  neglect  of  the  bottom  atreaa  (t^  ■  ■  0) ,  and  by 

the  praacrlptlon  of  a  purely  tonal  aurface  wind  atreaa.  At  the  rigid 
(and  aaaumed  vertical)  lateral  walla  of  the  bealn  a  tero-allp  boundary 
condition  u  »  v  ■  0  la  required,  from  which  the  tendency  of  the  water 
height  at  the  boundary  may  be  determined  from  the  continuity  equation, 
(21) ,  in  which  that  part  of  the  transport  divergence  involving  a  der¬ 
ivative  tangential  to  the  boundary  will  now  vaniah. 
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III.  THE  NUMERICAL  MODEL 

The  equations  (19)  to  (21)  are  approximated  on  a  finite-difference 
grid  scheme  in  the  surface  spherical  coordinates  X  and  ♦  as  sketched 
in  Fig.  1.  Of  the  many  possible  ways  to  assign  the  variables  on  the 
grid  points,  the  method  selected  here  involves  the  determination  of  the 
transport  u  and  v  at  points  for  which  (i  +  J)  is  odd,  with  the  water 
elevation  £  determined  at  the  interlocking  network  of  points  for  which 
(i  +  J)  is  even.  Here  i  and  j  are  integer  indices  for  the  longitude 
X  and  latitude  $  relative  to  the  origin,  such  that  X  ■  (i  -  l)AX/2, 

♦  ■  (J  -  1) A$/2 ,  with  i,j  -  1,2,...,  and  AX/2  and  A$/2  the  separations 
between  adjacent  grid  points.  Such  a  scheme  has  proven  useful  in  pre¬ 
vious  studies  (Cates,  1968),  and  is  particularly  convenient  with  the 
zero-slip  boundary  condition. 

Using  centered  differences  in  both  space  and  time,  and  neglecting 
the  bottom  stresses,  the  difference  equations  for  the  system  (19)  to 
(21)  are: 
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Fig.  1.  A  portion  of  the  spherical  finite-difference  grid  in  the 
reference  cm  A#  -  AX  -  2*,  shoving  ths  dsnslty  of  grid 
points  on  ths  ssas  seals  as  subsequent  figures.  Bara  the 
longitude  (X)  and  latitude  (♦)  of  a  grid  point  (relative 
to  the  origin  X  ■  0,  ♦  •  f  )  Is  given  by  ■  (1  -  l)AX/2 
and  i  •  f  ♦  (J  -  l)A+/2,  with  1  •  1,2,...,  61,  J  ■  1,2,..., 
49.  The  ?nset  shove  the  structure  of  the  Interlocking 
transport  (u,v)  and  elevation  (O  points. 
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Here,  for  example,  the  vertically  integrated  spaed  ■  u(lAX/2,  j At/2, 
nAt),  and  similarly  for  v°j  and  •  where  n  -  1,2,...  and  At  is  the 
tine  step  size.  At  the  initial  time  (corresponding  to  n  -  0) ,  a  single 
forward  time  step  ovsr  At/2  la  taken,  and  all  variables  referenced  at 
(n  -  1)  in  (22)  to  (24)  are  replaced  by  the  corresponding  variables  el 
n  (-  0).  In  the  nonlinear  tens  of  (22)  and  (23)  the  notation 

°?j  --IjV*  •  ^  (25) 

has  been  introduced  for  convenience,  since  the  water  depth  is  needed 
only  et  those  points  for  which  the  (vertically  integrated)  speeds  u 
and  v  ere  known.  The  frictional  tens  in  A  have  been  evaluated  from 
the  four  corner  or  diagonal  points  with  respect  to  ij  and  et  the  "back¬ 
ward"  ties  step  (n  -  1)  in  order  to  achieve  linear  computational  stabil¬ 
ity  with  the  leapfrog  differencing  scheme,  and  the  advectlvs  tense  in 
S  have  been  approximated  in  such  e  way  to  avoid  nonlinear  cospu rational 
instability  (Arekawa,  1966).  We  may  note  that  in  this  scheme  no  expli¬ 
cit  smoothing  or  interpolation  is  introduced. 

Since  (22)  and  (23)  ere  to  be  applied  only  et  points  for  which 
(1  +  J)  is  odd,  while  (24)  applies  only  to  points  for  which  (1  +  J)  la 
even,  the  necessary  variables  ere  automatically  provided  et  the  required 
grid  points.  A  scheme  in  which  all  three  variables  ere  determined  et 
each  grid  point  would  provide  solutions  no  more  accurate  than  those  from 
the  present  procedure  (for  the  same  AX,  A+,  and  At).  On  the  other  hand, 
e  scheme  in  which  u  is  determined  only  et  points  for  which  1  la  odd  and 
J  is  even,  end  v  determined  et  points  for  which  1  is  even  and  j  is  odd, 
would  be  e  more  economical  scheme  of  comparable  accuracy.  There  is  some 
tendency,  in  feet,  for  the  present  solutions  to  split  into  two  such  par¬ 
tial  solutions,  although  in  general  it  makes  little  difference  which  set 
is  selected. 
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In  all  of  tha  preaant  solutions  tha  ocaan  basin  is  rectangular  and 
of  unifora  depth,  but  in  order  to  accosnodate  later  studies  of  variable 
depth  (and  irregular  shape)  the  systea  (22)  to  (24)  has  not  been  rewrit¬ 
ten.  A  further  assuaptlen  in  the  present  work  (with  one  exception  not¬ 
ed  later)  ia  the  prescription  of  a  purely  sonal  wind  atress 


(  3*(J  ~  1)  \ 
X  *  —  \  2(J  -  1)  / 


-  -T  cos 


(24) 


-2 

where  T  ■  2  dynes  ca  and  where  j  ■  J  denotes  the  northern  wall  of  the 
baaln. 

The  cosrputatlonal  boundary  conditions  on  the  basin's  edges  aay  be 
written  as 


n  n  . 

**8  *  VB  "  0 


(27) 


for  all  n  and  for  all  boundaries  (where  B  denotes  boundary  points  fsr 
which  i  •  1,1  or  J  -  1,J),  and 
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for  n  ~  1,2,***.  Here  (28a)  to  (28d)  apply  to  the  western,  eastern, 
southern,  and  northern  boundaries,  respectively,  and  represent  the  ap¬ 
plication  of  the  continuity  equation,  (24),  at  the  boundary,  with  the 
centered  space  difference  replaced  by  a  one-sided  difference  between 
the  boundary  and  the  next  Interior  point,  together  with  the  condition 
(27).  At  n  ■  0  these  relations  are  replaced  by  the  corresponding  for¬ 
ward  time-difference  expressions. 

-2 

The  fixed  constants  have  been  taken  as  a  -  6371  ka,  g  -  980  cm  sec  , 
-5  -1  -3 

n  -  7.292  *  10  sec  ,  and  pQ  ■  1  g  c«  ,  while  the  reaalnlng  constante, 
fixed  for  any  particular  integration,  are  varied  in  a  series  of  compara¬ 
tive  solutions.  A  set  of  typical  values  (which  specify  a  reference 
8  2  —1 

case)  are  A  ■  10  cm  sec  ,  S  ■  4/3,  +o  (southern  boundary)  “  18  N, 

AX  ■  A+  -  2  deg.  At  ■*  10  min,  and  the  basin's  zonal  and  meridional  widths 
(I  -  l)AX/2  and  (J  -  l)A+/2  given  with  I  ■  61,  J  •  49.  The  integra¬ 
tions  have  proved  stable  when  carried  out  for  as  long  as  60  days  of 
oceanic  (simulated)  time  starting  from  an  initially  undisturbed  ocean 
(u°  ■  v°  ■  C°  “  0),  and  require  approximately  1  min  of  machine  time  per 
day  of  simulated  time  on  the  CDC  6600  computer. 
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7V.  C0MPA1ATIVE  NUMERICAL  SOLUTIONS— NUMERICAL 
AND  GEOMETRICAL  EFFECTS 


RESOLUTION  TESTS 

In  a  program  of  comparative  numerical  Integration,  one  of  the 

first  questions  to  be  exaained  should  be  the  level  of  truncation  error 

in  the  solutions.  This  is  particularly  important  in  the  present  prob- 

lea  in  view  of  the  narrowness  of  the  western  boundary  currents  and 

the  characteristic  presence  of  short  planetary  transient  waves.  We 

thus  first  examine  the  solutions'  behavior  es  the  resolution  of  the 

finite-difference  grid  is  systematically  changed,  while  all  other  para* 

meters  are  held  constant.  For  this  purpose  it  is  sufflcelnt  to  examine 

only  the  solutions  for  meridional  transport,  generated  in  the  basin  of 

fixed  dimensions  and  subject  to  the  steady  wind  stress,  (26) .  This 

configuration  will  constitute  a  reference  case,  for  which  we  select 

AX  (I  -  l)/2  ■  60  deg  end  A+(J  -  l)/2  -  48  deg  as  the  zonal  and  marl- 

dlonal  dimensions  of  the  basin,  h  -  400  m,  S  •  4/3,  T  -  2  dynes  cm'  , 

8  2  -1 

4  ■  18  N,  and  A  ■  2  *  10  cm  sec  Under  theae  conditions  the  solu- 
tions  of  (22)  to  (24)  for  the  meridional  transport  v  at  30  days  are 
shown  in  Fig.  2  for  A+("AA)  ■  3*  and  2*,  and  each  is  compared  with 
the  higher-resolution  case  A4(*AA)  ■  3/2*. 

Although  the  general  character  of  the  (meridional)  circulation  is 
similar  in  the  3*  and  3/2*  cases  (Fig  2a),  the  3*  solution  shows  a 
systematic  eastward  displacement  (relative  to  the  3/2*  case)  of  the 
western  countercurrents  and  the  mid-ocean  transients.  Since  these  fea¬ 
tures  evolve  as  westward-propagating  lossby  waves,  we  may  Interpret  this 
difference  ae  a  phase  truncation  error.  In  the  western  boundary  currents 
themselves,  the  3*  solution  has  overestimated  the  maximum  of  the  north¬ 
ward  current  and  similarly  underestimated  the  southward  current,  as 
stmmarlted  in  Table  1,  and  in  addition  has  positioned  the  flow  too  close 
to  the  boundary.  In  the  southeastern  portion  of  the  basin  a  spurious 
southward  "boundary"  current  le  also  present  in  the  3*  solution,  and 
persists  throughout  the  calculation.  This  feature  is  apparently  rela¬ 
ted  to  the  solution's  overestimate  of  the  northward  flow  in  the  western 


Fig.  2.  The  solutloaa  for  mtUIomI  trioi^ort  v  U  Sr  (100  ka) 
la  a  spherical  basin  at  30  days  for  M(*  AX)  •  3*  (a) 

•ad  A#(«  AX)  >2*  (b) ,  coopered  with  cha  aolutloo  far 
A4(«  AX)  •  3/2*  shown  by  dashed  lines  ta  both  (a)  sad  (b) 
Bara  A  -  2  *  10*  cm2  sac"1,  h  •  400  a;  the  baaia  is  40 
vida  and  sat ends  fraa  It  H  ta  44  H.  The  (vertically -aver 
aged)  currant  v/h  any  be  faaad  la  ca  aac”l  by  aaltlplylag 
the  laoliaaa  by  tba  factor  3/2. 
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perc  et  the  Natktrv  gyre  (m«  Table  1)#  m4  coetrlbetes  to  tho  rela¬ 
tively  high  total  ImU  klootlc  raergy  of  i.l  ■  10  erge,  coopered 
vlth  4.0  ■  1023  ergs  fer  the  helved  0/2*)  grid.  le  overall  perfora- 
eace  the  3*  grid  therefore  dees  aec  provide  saff Iciest  resolutloe. 

The  coaperlsoe  of  the  2*  ead  3/2*  soletloae  la  Fig.  2b,  oa  the 
other  head,  sheas  soasohat  better  cerre speeds oce .  There  Is  little  aya- 
teaetlc  phaee  d If foresee  evldeat,  aad  the  aagaltedee  of  the  stroager 
trsas ports  la  the  asst  ere  la  each  closer  agreeaeat  thaa  la  Fig  2a  (aae 
alas  Table  1).  The  tmacatlee-ladoced  aa seal  ass  cs treat  aloag  the  Mat¬ 
ed  shore  Is  also  she eat  vlth  the  2*  (aad  3/2*)  grid,  aad  the  baa la' e 
total  klaatlc  eaergy  of  4.1  ■  1023  ergs  la  la  good  agreeaeat  vlth  that 
for  the  3/2*  caee.  Oa  this  evldeece  ae  aay  coaclede  that  the  2*  grid 
gleee  adages te  reoolvtloa  through oet  the  haala  aad  repreeoate  aa 
eealc  eelectloa  a a  e  re foresee  reeoletloa  la  farther  ceaparetlve 
gretloaa.  Dee  of  e  flaer  grid  la  the  vast  era  part  af  the  baela, 
ever,  avoid  revolt  la  e  slightly  better  reeelatlea  of  the  lateas 
ary  c arrests  (aa  recegalsed  earlier  by  iryea  (1H3)) ,  altheogh  tl 
of  the  tmacatlea  errwr  has  bees  reaevad  la  paeelag  frea  the  3* 

2*  grid. 

Me  aay  see  the  data  ef  Table  1  te  eat last a  the  obsolete  else  af 
the  aolatleaa '  tmacatlea  error  ea  the  as  empties  that  It  varies  as  the 
oguare  ef  the  spatial  aesh  alas.  I  lace  the  aa  visas  carreata  la  the  seat 
gees rally  de  aet  occur  et  celacldeat  palate  ef  the  several  gride,  the 
total  or  sat  traaapert  la  the  ousters  beuadary  ce treat  eervee  aa  a  cos¬ 
parable  sees  ere  of  tmacatlea  error.  The  aet  sort  hoard  beuadary  curreat 
traaapert  a  la  the  sou  there  aad  earths  ru  gyres  frea  Table  1  am  thus 
20  tv  (Iv  •  10*  a3  aec’1),  13  tv,  aad  14  9v  fer  the  3*,  2*,  aad  3/2* 
grid  solo t leas.  If  these  traaoperte  ceatala  a  tmacatlea  error  p ro¬ 
per  tloaal  te  (Ag) 2  superposed  epee  the  true  or  trvacet loe-f ree  value, 
thee  the  ratle  (Hj  -  HjHWj  -  H,^)"1,  fer  esaaple,  should  egual  20/27, 
ahem  le  the  aet  aerthverd  beuadary  curreat  traaapert  la  the  A#  grid 
aelutlea.  Vlth  the  above  values  this  ratle  la  3/A,  vtilch  aay  be  takes 
aa  reaeoaable  ceaf  lraatlea  ef  the  tmacatlea  hype  thesis.  If  M  deaotea 
the  tmacat lea-free  selotlea,  thee  a  perceataga  tmacatlea  error  aay  be 
avpreaaed  by  (N^  •  M^M*1,  fer  esaaple,  glvlag  e  23  pe reset  error  for 


» 


•la¬ 


th*  2*  solution  for  th*  western  boundary  current  (when  th*  value  M  ■ 

o 

12  Sv  a*  determined  by  the  3*  end  3/2*  grids  Is  used) .  This  Is  e  max¬ 
imum  error  for  th*  2*  grid,  as  the  transport  of  the  longer  reflected 
transient  waves  In  the  Interior  (sea  Fig.  2)  ere  determined  with  some¬ 
what  greater  accuracy. 


SPHERICAL  VERSUS  6 -FLAKE  COORDINATES 


In  view  of  th*  widespread  use  of  th*  8-plen*  approximation  In 
oceanic  modeling.  It  Is  of  Interest  to  compare  such  solutions  with 
those  from  th*  present  spherical  formulation.  Th*  6-plane  dynevl 
equations  may  be  obtained  from  the  system  (19)  to  (21)  by  replacing 
cos  4  by  unity  throughout,  neglecting  the  terms  In  ten  #,  and  writing 
dx  “  ad A,  dy  “  ad4.  The  resulting  equations  are: 


Ju 
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where  u  and  v  continue  to  ba  the  vertically  Integrated  volume  transports 
(u)  and  (v),  and  V*  le  the  Laplaclan  In  the  rectangular  coordinates  x  and 
y.  Hare  th*  Corloll*  parameter  le  now  given  by  f  -  f  ♦  6y,  with  fQ  and 
6  taken  as  constant*.  For  a  uniform-depth  basin  (and  with  S  •  1)  thi* 
eystem  1*  the  earns  as  that  used  In  a  previous  etudy  (Cat**,  1968)  when 
tha  meridional  wind  etree*  component  and  both  component*  of  the  bottom 
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•treas  arc  neglected.  Thla  B-plana  system  may  alao  ba  found  by  a  for¬ 
mal  tr  ana  format  Ion  of  tha  apbarlcal  equations,  aa  shown,  for  tuapla, 
by  Varonla  (1963),  and  la  valid  for  tha  quasl-herlsontal  mot Iona  In  a 
•hallow  ocaan.  Tha  dlffaranca  equations  for  tha  Bflana  caaa  may  ba 
raadlly  writ tan  by  appropriate  alapllf lcatlon  of  tha  apharlcal  aqua¬ 
tions  ,  (22)  to  (24),  and  tha  (-boundary  condition  (28). 

In  tha  praaant  application  va  continue  to  Ignore  tha  bottom  atraaa, 

and  aalact  tha  constant*  A  ■  10®  cm*  a»c’\  h  ■  400  m,  S  •  4/3,  I  ■  61, 

and  J  •  49  for  both  tha  6*?lana  and  apharlcal  casaa,  with  tha  wind 

atraaa  glvan  by  (26).  For  a  rafaranca  aolutlon  an  tha  aphara  wa  aalact 

AX  •  44  •  2*  and  4  -  18  H;  for  tha  comparative  solutlona  on  tha  8* 

plana  wa  select  Ax  •  Ay  •  222  ka,  f  ■  4.307  *  10  3  sac*  ,  and  8  • 

—13  —1  —1  ® 

1.701  *  10  cm  aac  ,  tha  latter  corresponding  to  the  baaln'a  aver¬ 
age  latitude,  42  R.  Theaa  comparative  aolutlona  are  shown  In  Flga,  3 
and  4,  and  aarva  at  tha  same  time  to  document  tha  time -depandent  behav¬ 
ior  In  tha  apharlcal  caaa. 

Tha  evolution  of  tha  meridional  tr ana port,  v,  shown  In  Fig.  3  (solid 
lines)  consists  of  a  aarlas  of  transient  Roaaby  waves  reflected  from 
tha  western  shore,  with  these  waves  undergoing  progressiva  slowing, 
shortening,  and  amplitude  decay  aa  they  marge  with  tha  quaal-eteady 
boundary  currents.  These  transients  are  particularly  prominent  In  tha 
southern  gyre,  where  they  display  a  characteristic  wavelength  of  soma 
13  to  20  dag  (longltuda),  and  move  westward  at  speeds  of  approximately 
1.2  dag  longltuda  day'*  (about  1.3  a  sac**  at  tha  gyre*a  center,  34  R). 
While  tha  western  boundary  current  la  In  apparent  equilibrium  before 
tha  21st  day,  tha  countercurrent  offshore  In  the  southern  gyre  responds 
to  tha  periodic  arrival  and  absorption  of  tha  transients;  the  south¬ 
ward  transport  In  tha  west-central  basin  at  21  days,  for  axaaple,  aay 
ba  aeon  merging  with  the  aoutharn  gyra's  countercurrent  at  30  days, 
uhlla  a  succeeding  transient  now  occupies  the  former's  offshore  posi¬ 
tion.  It  alao  appears  that  by  48  days  thasa  transients  have  undergone 
reflection  from  tha  southern  boundary  of  tha  basin  (Cates,  1969).  In 
tha  northern  gyre  tha  (northward)  countercurrent  Is  In  position  after 
•bout  30  days,  and  shows  little  subsequent  change.  The  transient  Ross- 
by  waves  In  this  gyrs  ars  also  seen  to  move  slowly  to  tha  west  (averaging 


The  solutions  tor  meridional  transport  v  in  Sv  (100  ka) 
in  a  spherical  baaln  with  AX  •  A#  -  2*  (solid  lines)  and  in 
a  6-plane  with  Ax  -  Ay  •  222  ka  (dashed  lines) .  Here  A  - 
10®  ca*  sec~l,  h  -  400  a  and  the  basin  extends  froa  18  N  to 
66  N.  The  saae  isoline  spacing  (uneven)  has  been  used  for 
both  the  spherical  and  6-plane  solutions!  with  only  the 
sero  isoline  labeled  (in  parentheses)  in  the  latter  case. 
The  6-plane  solutions  have  been  spaced  eastward  froa  the 
western  boundary  according  to  sec  4  to  approxiaate  the 
geoaatry  of  the  spherical  solutions.  The  aean  current  v/h 
in  ca  sec~l  is  given  by  aultlplylng  the  transport  lsolines 
by  factor  5/2. 
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about  0.4  dag  longltuda  day  *),  and  ara  less  veil  developed  than  their 
counterparts  In  the  southern  gyre.  These  sase  general  features  are 
also  reflected  In  the  solutions  for  the  tonal  transport,  u,  and  surface 
elevation,  C»  shown  In  Fig.  4  at  30  days,  although  the  role  of  the 
transients  Is  less  conspicuous  In  these  variables. 

Of  particular  lnterast  are  the  comparative  solutions  for  the  6- 
plane  also  shown  In  Figs.  3  and  4  (dashed  lines).  Here  the  northward 
finite-difference  mesh  size  Is  the  same  In  the  spherical  and  B-plane 
solutions  (2*),  while  the  eestvsrd  mesh  size  Is  different  at  each  lat- 
tltude  over  the  basin  (2*  longitude  compared  with  222  km) .  Starting 
with  points  coincident  with  the  spherical  grid  on  the  western  boundary, 
the  solutions  on  the  B-plane  have  here  been  plotted  with  the  points  at 
each  latitude,  4,  displaced  eastward  according  to  2*  sec  4,  In  order  to 
remove  the  zonal  scale  distortion  of  the  spherical  geometry.  The  east¬ 
ern  portion  of  the  B-plane  grid  Is  therefore  not  seen  In  these  figures, 
and  the  two  solutions  should  not  be  compared  In  the  eastern  third  or 
so  of  the  (60 -dag)  grid  shown. 

In  the  western  and  central  portions  of  the  basin,  we  see  a  close 
correspondence  In  the  general  patterns  of  the  two  solutions,  although 
there  are  two  significant  differences.  In  the  southern  gyre  the  trans¬ 
ients  of  the  B-plane  solution  lie  systematically  to  the  east  of  their 
spherical  counterparts.  This  displacement,  moreover,  Is  cumulative  with 
time,  and  Indicates  that  the  westward  phase  speed  of  the  B-plaue  trans¬ 
ients  Is  approximately  90  percent  of  that  In  the  spherical  solutions. 
This  effect  Is  evidently  not  due  to  the  differences  In  zonal  truncation 
error  In  the  two  solutions,  as  the  phase  differences  are  smeller  In  the 
northern  gyre  where  the  solutions'  zonal  mesh  discrepancy  Is  the  great¬ 
est.  A  second  systematic  difference  Is  a  northward  displacement  of  the 
B-plane  solutions  of  about  2  deg  relative  to  the  spherical  case.  This 
effect  Is  present  over  most  of  the  basin  and  persists  throughout  the 

calculation,  and  may  be  related  to  the  choice  of  the  constants  f  and  6. 

o 

Even  a  cursory  examination  of  Figs.  3  and  4,  however,  shows  that 
the  greatest  difference  between  the  spherical  and  B-plane  solutions  is 
the  somewhat  larger  amplitude  of  the  transports  (and  water  elevation) 

In  the  6-plane  case:  In  general  this  excess  amounts  to  about  SO  percent 


Tha  solutions  for  tonal  transport  u  in  Sv  (100  ka)  (a) 
and  for  vatar  slsratlon  (  In  ea  (b),  accompanying  tha 
solutions  of  Fig.  3  at  30  days.  Tbs  dashsd  Unas  ars  tbs 
solutions  on  tbs  6 -plana  (sss  Fig.  3). 
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of  cha  spherical  case's  amplitudes.  More  detailed  data  on  theae  dif¬ 
ferences  are  given  in  Table  2  for  the  Important  western  boundary  cur¬ 
rent,  which  we  note  to  be  generally  stronger  and  alightly  narrower  in 
the  6-plane  case.  Since  the  boundary  conditions,  basin  depth,  lateral 
viscosity,  marldlonal  grid  site,  and  the  lmpoaed  wind  atresa  are  the 
same  In  both  solutions,  the  source  of  thia  difference  la  to  be  found 
In  the  geometry  of  the  6-plane  approximation  Itself;  namely,  the  pro¬ 
gressive  narrowing  with  latitude  of  the  basin's  zonal  width  in  the 
spherical  formulation  with  meridional  boundaries.  Although  the  merid¬ 
ional  dimensions  of  the  two  basins  are  the  same,  the  sonal  dimensions 
differ  by  more  then  a  factor  of  two  at  the  northern  boundary.  With  a 
Sverdrup  equlllbrlun  over  the  Interior  of  the  basin,  the  Integral  of 
the  wind-stress  curl  Is  proportional  to  the  total  interior  meridional 
flow,  and  hence,  through  continuity,  is  a  measure  of  the  net  flow  in 
the  western  boundary  currants.  In  the  present  cases,  the  ratio  of  the 
area  Integral  of  the  curl  of  the  wind  streae  for  the  epherlcal  and  6- 
plane  basins  Is  0.42,  and  the  ratio  of  the  net  northward  traneport  In 
the  western  boundary  currents  In  the  two  cases  is  0.48  from  the  data 
of  Table  2.  This  agreement  would  be  even  closer  were  the  western 
countercurrents'  transport  also  considered.  The  magnitude  dlfferencea 
In  the  eolutlona  of  Figs.  3  and  4  are  thus  believed  to  be  due  to  this 
areal  effect,  and  close  qusntltatlve  agreement  could  be  expected  if 
the  eastern  boundary  In  the  6-plane  basin  were  moved  progressively 
westward  with  latitude  In  order  to  approximate  the  sonal  dimension  of 
the  spherical  basin. 

In  theoretical  studies  of  the  adequacy  of  the  6-plane  approxima¬ 
tion  In  a  nonvlscous  model,  Longue t-Hlgglns  (1964,  1965)  has  shown  that 
the  periods  of  at  least  the  lower-order  modes  (of  free  oscillations)  In 
simple  enclosed  basins  of  modest  else  should  agree  to  within  about  10 
percent  In  the  spherical  and  6 -plane  solutions.  This  conclusion  is 
supported  by  the  transient  behavior  discussed  above  In  connection  with 
Fig.  3.  From  analytical  solutions  for  a  linear,  nondlvergent  model 
with  bottom  friction,  Takano  (1966)  has  also  shown  that  the  principal 
difference  between  (steady)  spherical  and  6-plane  solutions  Is  the  in¬ 
creased  boundary  current  transport  of  the  latter  at  higher  latitudes, 
due  to  the  Increased  basin  width  of  the  6-plane  formulation.  This  is 


confined  by  the  present  solutions,  which  nay  therefore  be  taken  as  a 
demonstration  of  the  adequacy  of  the  8 -plane  approximation  for  the 
time-dependent  motions  In  the  lower  and  middle  latitudes.  In  the  num¬ 
erical  simulation  of  realistic  ocsanlc  motions,  however,  there  seems 
to  be  no  reason  not  to  employ  the  more  general  spherical  system,  end 
It  will  be  used  In  the  furthsr  comparative  solutions  discussed  below. 

EFFECTS  OF  THE  BOUNDARY  AND  BASIN  SIZE 

In  a  wind-driven  oceanic  model  the  location  of  the  regions  of 
zonal  Inflow  and  outflow  to  the  western  boundary  current  is  determined 
both  by  the  spatial  variation  of  the  Imposed  wind  stress  and  by  the 
location  of  the  zonal  boundaries  of  the  basin.  In  an  effort  to  docu¬ 
ment  this  latter,  geometrical  effect  as  distinguished  from  the  former, 
dynamical  one,  a  comparative  lntegretlon  was  made  in  which  the  southern 
boundary  was  moved  farther  south,  as  well  as  a  second  Integration  in 
which  the  basin's  zonal  width  was  reduced.  The  same  zonal  wind-stress 
pattern  was  used  In  each  case,  with  the  solutions  of  a  (spherical)  case 
of  Intermediate  basin  size  serving  as  a  reference. 

With  the  constants  A  «  2  *  10®  cm2  sec  \  AX  •  A#  ■  2*,  h  •  400  m 
and  S  ■  4/3,  the  solutions  for  the  transports  at  30  days  are  shown  in 
Fig.  5  for  the  comparative  b sains  of  48  deg  and  64  deg  latitudinal  ex¬ 
tent.  The  wind  stress  la  given  by  (26)  for  the  48-deg  case,  and  by 
r*  •  -T  sin  [2w(J  -  1) (J  -  l)"1],  t*  ■  0  for  the  64-deg  case  (with  J  •  1 

now  at  2  N)  to  ensure  coincidence  of  the  stress  In  the  region  common  to 

„2 

both  basins.  (In  all  cases  T  ■  2  dynes  cm  .)  The  Imposition  of  the 
southern  boundary  at  18  N  results  In  a  dscrsase  of  the  westward  and 
northward  transports  near  the  southern  and  western  boundaries,  respec¬ 
tively,  as  shown  by  the  dashed  lines  in  Fig.  5.  Beyond  this  boundary 
zone  of  some  2  to  4  deg  width,  the  presence  of  the  southern  boundary  at 
18  N  produces  an  Increased  westward  transport  relative  to  the  2  N  case 
In  a  broad  zonal  band  extending  approximately  half  way  across  the  basin. 

The  differences  In  the  meridional  transport,  also  shown  In  Fig.  5. 
are  similarly  confined  to  the  southern  and  western  regions  of  the  basin. 
Imposition  of  the  southern  boundary  at  18  N  (rather  than  at  2  N)  has 
generally  limited  the  southward  transport  In  the  Interior  of  the  basin, 
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Tha  aolutlona  for  tonal  traaoport  (•)  and  narldlonal 
tr  ana  port  (b)  In  8r(100  ka)~l  at  30  daya  vlth  tha 
•ou than  boundary  raaovad  to  2  If.  Bara  AX  -  A#  -  2#, 

A  ■  2  x  10®  cn*  aac"1,  and  h  ■  400  a.  Tha  dlffarancaa 
with  raapact  to  eorraapondlng  aolutlona  vlth  tha  aoutharn 
boundary  at  16  M  ara  ahovn  by  tha  daahad  llaaa  (2  If  eaaa 
■lnua  16  N  caaa) .  Tha  (vartically-avaragad)  curran ta 
u/h  and  r/h  in  cn  aac”*  ara  glvan  by  aultlplylng  tha 
tranaport  laollnaa  by  tha  factor  3/2. 


nd  limited  the  northward  transport  la  nearly  all  of  the  wee  tern  bound¬ 
ary  currant  eye  tea.  The  largest  boundary  af  facta  oa  both  the  tonal  nd 
meridional  trasaport  are  thua  found  where  the  flow  ia  the  noat  lntanaa , 
which  iodteataa  that  the  affect  of  a  particular  boundary  ia  wary  much 
dependant  upon  ita  location  with  raapect  to  circulation  inferred  from 
the  driving  atraoa  pattern.  The  II  N  boundary  in  Fig*  3  ia  at  a  lati¬ 
tude  of  email  meridional  tr  ana  port  but  large  tonal  traarport  in  tha 
2  H  solution;  ware  tha  boundary  at  34  N,  for  exanple,  laao  affect  on 
tha  tonal  flow  but  more  affect  on  tha  meridional  flow  could  be  expected 
In  tha  aouthern  interior  of  tha  II  N  baa  in  tha  boundary-induced  trano- 
port  chengee  are  of  tha  order  of  23  to  30  percent  of  tha  eolwtiona  with 
a  wore  diatent  aouthern  boundary,  and  th la  affect  ahould  be  kept  In 
mind  whan  interpreting  near-boundary  clrcwlatloae .  A  eimller  aenoitlw- 
ity  to  relocation  of  tha  northern  boundary  could  elao  be  expected  In 
the  preeent  caaae. 

In  a  subsequent  experiment  the  location  of  the  oaetern  boundary 
waa  changed  from  that  in  Fig.  3  to  a  poeitlon  in  the  center  of  thie 
baa in.  The  tranaport  change a  produced  at  30  daya  by  thla  alteration 
are  ahown  in  Fig.  6,  again  referenced  to  the  larger  (10-deg  long) 
baa  in' a  aolutiona.  The  contribution  of  the  aaatara  half -baa  in  la  aeon 
to  be  reeponeible  for  approximately  half  of  the  naxlman  trana porta  In 
the  weatern  ocean.  In  particular,  the  meridional  wee  tarn  boundary  cur¬ 
rent*  are  preeent  on  the  a  ana  a  cal  a  in  both  solutions,  and  are  wary 
nearly  halved  in  a  peed  by  the  halving  of  the  baola'e  aaat-weet  dlnan 
a  loo ,  while  the  maxinmn  tonal  tranaport  ie  reduced  by  approximately  10 
percent.  Since  the  a  ana  tonal  wind-etreoa  dietrlbwtioa  wee  impound  In 
both  caaea,  theee  tranaport  chaogee  ere  the  direct  reault  of  the  baain' 
changed  width.  Tha  interior  meridional  flow  la  little  affected  by  tho 
width  chaagee,  and  the  net  meridional  tranaport  accwmelatee  with  dla- 
tance  from  the  eaatern  boundary  aa  ahown  by  Sverdrup  (1947).  Tho  la¬ 
tent  lty  (but  net  the  ecale)  of  the  oaetern  boundary  current  ie  then 
determined  by  the  continuity  requirement  to  balance  thla  flow.  Tho 
total  kinetic  energy  in  the  30-  and  10-deg  longitude  baalaa  at  30  days 
la  1.0  "  10**  «rga  and  3.1  ■  10**  urge,  reepectlvaly,  and  la  e  further 
reflection  ef  the  approximate  halving  of  the  larger  trana  porta.  Thla 
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Tha  oolatloat  for  mmI  transport  (a)  aa4  aorl41oaal  trans¬ 
port  (b)  la  Sr  (100  la)’1  at  30  days,  with  ths  daehad  llaaa 
ladl catlap  tba  dlfforaara  la  transport  with  the  aaotara 
Soiiary  at  30  dap  leapIfaOe  rather  than  at  SO  dan  laofltuda 
(30  dap  caaa  alaas  SO  dap  caaa) .  Bara  A  •  2  *  10*  cr  sac"’ 
h  •  S00  a,  aai  41  •  AS  •  2*  as  la  Ftp.  3.  Tba  ararapa  enr- 
raata  la  ca  sac"1  ara  plvea  bp  aaltlplplop  tba  traaapart 
lao llaaa  bp  tba  factor  3/2. 
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hh  affect  wii  prrrleuljr  14ntlflH  m  the  likely  causa  sf  ths  eya- 
tanotlc  nognltuda  differences  between  Cha  otherwise  ceapareble  solutions 
os  tba  sphere  sad  6 -plane  Is  Figs.  3  and  4. 

Thu  boundary  Indue ad  circulation  changes  discus sad  above  should 
sot  ba  rstsrdad  aa  "errors"  la  say  sasaa,  but  ratbar  aa  Indications 
of  tha  sensitivity  of  ths  solutions  to  ths  actual  placsasat  of  tha 
basis's  lstsrsl  boundaries.  Aa  applied  to  ths  world's  oceans ,  these 
results  would  suggest  that  tha  vertically- 1st agratad  transport  of  tha 
western  boundary  currents  la  ths  Pacific  ocera,  for  aiaayla,  should 
bo  approx lawtsly  double  those  In  ths  Atlantic ,  la  else  of  ths  two  acaaaa 
approx  lasts  2  to  1  tonal  d  leans  Ion  ratio  and  their  generally  ceuparabla 
stress  fields.  Tha  fact  that  tbs  Kwroehio's  trees port  is  gone rally 
lass  than  that  of  the  Oulf  St  rasa  any  be  due  to  ths  preaeece  of  a  large 
borocllale  coup  onset  In  the  letter's  transport,  or  to  ether  factors  not 
considered  In  tha  present  wlad -driven  homogeneous  nodal. 
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V.  CCHPAIATIVE  WOKERICAL  SOLUTIONS— PHYSICAL  EFFECTS 

FKICTIOWAL  gmcrs 

Although  the  wind-driven  circulation  In  the  Interior  of  an  ocean 
baa  In  la  generally  characterised  by  a  Sverdrup  equilibrium  (approxi¬ 
mate  balance  between  the  local  wlnd-etraaa  curl  and  the  planetary  vor- 
tlclty  (0)  tendency),  the  western  boundary  currants  are  known  to  be 
critically  Influenced  by  frictional  forces •  While  the  boundary  cur¬ 
rants  may  not  owe  their  existence  In  a  fundamental  sense  to  the  fric¬ 
tion  (when  we  view  such  currents  as  tha  result  of  the  trapping  of  west¬ 
ward  propagating  Rossby  waves) ,  their  width  end  Intensity  are  prin¬ 
cipally  datarmload  by  viscous  offsets  If  we  assume  that  tha  nonlinear 
Inertial  affects  are  not  dominant. 

Since  the  work  of  Hunk  (19S0)  a  lataral  eddy  viscosity  A  of  the 
7  12-1 

order  10  to  10  cm  sec  In  wind-driven  models  has  provided  reason¬ 
able  values  for  the  width  and  transport  of  tha  western  boundary  cur¬ 
rent,  although  actual  boundary  currents  sppesr  to  be  systematically 
narrower  and  more  Intense.  The  variation  of  such  solutions  to  changes 
In  tha  lateral  viscosity  has  been  documented  by  Bryan  (1963)  on  the 
0 -plane :  a  relatively  large  value  of  A  will  produce  a  smoothed  and 
"sluggish"  circulation,  while  a  relatively  small  value  will  produce 
an  Intense  and  narrow  boundary  current  which  may  become  unstable.  Pre¬ 
sumably  a  similar  behavior  would  occur  in  the  present  (spherical)  for¬ 
mulation,  and  aome  data  are  given  In  Table  3;  here  solutions  have  not 

7  2-1 

been  accomplished  for  values  of  A  leas  than  6  *  10  cm  sec  with  the 
present  2 -deg  resolution.  These  data,  as  well  ss  a  comparison  of  Figs. 
2  and  3,  show  that  the  Intensity  of  the  meridional  current  at  30  days 
is  multiplied  by  a  factor  of  approximately  0.7  by  the  doubling  of  A 
from  10*  cm*  sec”1  to  2  *  10*  cm*  sec*1.  A  closer  inspection  shows 
that  the  western  boundary  current  in  both  the  northern  and  southern 
gyres  le  approximately  10  percent  wider  In  the  cass  of  the  larger  A. 
These  results  may  be  compared  with  the  factors  2  11  and  2^^  to  be 
applied  to  the  intensity  and  width,  respectively,  accompanying  a  doub¬ 
ling  of  A  according  to  the  theory  of  a  (linear)  viscous  boundary  layer 


emcTS  or  fkictiom  on  toe  wunut  soukdaiy  cum  ext* 
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(Hunk,  1950),  whereas  the  current's  totsl  transport  is  Independent 
of  A.  There  Is  elso  s  see 11  (2-deg  longitude)  but  systeaatic  eastward 
displacement  of  the  transients  in  the  centrsl  ocean  for  the  case  of  the 
larger  A,  as  well  as  s  slight  amplitude  reduction.  In  this  respect  the 
effect  of  increasing  the  lateral  viscosity  is  similar  to  that  produced 
by  an  Increase  in  the  grid-mesh  size  noted  earlier  in  connection  with 

rig.  2. 

The  viscous  dissipation  In  a  wind-driven  homogeneous  model  may 

also  be  formulatad  in  terms  of  an  Ekman  (or  bottom)  friction  in  the 

manner  of  StoMel  (1948).  Some  evidence  on  the  response  of  the  neer- 

equlllbrlum  circulation  of  such  an  ocean  to  variations  of  the  bottom 

friction  coefficient  has  been  given  by  Veronls  (1966b,  1966c),  but 

only  for  systems  In  which  the  nonlinear  inertial  terns  were  relatively 

mors  Important.  Since  this  model  provides  for  no  countercurrents  and 

characteristically  places  the  maximum  meridional  current  at  the  wsatern 

boundary  ltaelf,  integrations  have  not  been  made  with  the  replacement 

of  lateral  friction  by  bottom  friction  in  the  present  formulation. 

Bather,  the  sensitivity  of  the  prassnt  model  with  lateral  viscosity 

to  the  addition  of  a  bottom  frictional  dissipation  has  been  considered. 

The  results  of  such  a  comparative  Integration  are  shown  In  Fig.  7  for 

A  ■  10®  cm"'  sec  *  end  with  the  coefficient  k  -  10~®  sec  *  in  the  added 

bottom  frictional  terms  -ku  and  -kv  on  the  right-hand  sides  of  (19) 

and  (20),  respectively.  (These  terns  were  evaluated  at  the  backward 

time  level  n  -  1  In  the  finite-difference  equations  (22)  and  (23)  to 

promote  computational  stability.)  Also  shown  as  dashed  lines  in  Fig. 

7  ere  the  comparative  or  reference  solutions  for  the  case  of  a  lateral 
8  2  -1 

viscosity  A  ■  10  cm  sec  alone. 

The  maximum  zonal  flow  in  the  western  part  of  the  baaln  Is  seen 
to  be  reduced  by  the  bottom  friction  to  about  2/3  of  Its  value  with  lat¬ 
eral  friction  alone,  although  the  location  and  shape  of  the  circulation 
patterns  are  similar.  In  the  remainder  of  the  baeln  the  zonal  flow 
(which  Is  very  nearly  geostrophlc)  Is  affected  very  little  by  the  Intro¬ 
duction  of  bottom  friction.  This  indicates  that  the  north/south  pro¬ 
file  of  the  free  water  surface  le  maintained  primarily  by  the  meridio¬ 
nal  transport  In  the  surface  Ekman  layer  Induced  by  the  zonal  wind 


Tha  solutions  for  soosl  transport  (a)  and  aarldlonal  trans¬ 
port  (b)  In  Sv  (100  kn)"*  at  30  days  with  a  bo t ton  fric¬ 
tion  coefficient  k  ■  10"®  sac”*  and  lateral  viscosity 
A  -  10®  cn2  sac"*  (solid  lines) .  The  corresponding  solu¬ 
tions  with  lateral  viscosity  alona  are  shown  by  dashed 
lines.  The  nean  currants  u/h  and  v/h  In  cn  sac”*  era  given 
by  Multiplying  the  transport  Iso  Unas  by  tha  factor  5/2 
for  tha  present  400  a  basin. 
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stress  distribution.  The  maximum  meridional  current  in  Fig.  7  near 
the  western  boundary  is  also  reduced  to  approxlautely  two-thirds  its 
value  with  lateral  viscosity  alone,  while  the  current's  width  is  in¬ 
creased  by  the  introduction  of  bottom  friction.  From  these  solutions 
(and  the  dsta  of  Table  3)  we  deduce  that  the  value  k  -  10~^  sec  *  has 

approximately  the  same  effect  on  the  maximum  meridional  current  as  does 

8  2  -1 

e  value  of  A  of  about  10  cm  sec  Bottom  friction,  however,  reduces 
the  total  transport  of  the  western  boundary  current,  in  contrast  to 
the  transport's  Independence  of  the  lateral  viscosity  (with  no  bottom 
friction). 

This  dynamical  equivalence  applies  only  to  the  quasl-statlonary 
boundary  currents,  however.  We  note  from  Fig.  7  that  the  transients 
prominent  in  the  west-central  ocean  with  only  lateral  viscosity  have 
been  almost  completely  damped  out  by  the  introduction  of  bottom  fric¬ 
tion.  This  illustrates  the  sensitivity  of  transient  Rossby  waves  to. 

this  form  of  dissipation.  If  we  Interpret  the  magnitude  of  the  lat- 

2 

eral  viscous  term  AV  v  for  these  transients  in  terms  of  a  charscteris- 

2  -1 

tic  decay  time  t^  ■  v^tAV  v)  ,  where  v^  is  a  typical  transient's  am¬ 
plitude,  we  find  t.  ■  46  days  for  the  typical  wave  of  1600  km  length 

g  2  1  — 1 

when  A  ■  10  cm  sec  Since  the  period  of  these  waves  (about  14  days; 
see  Fig.  3)  is  somewhat  less  then  this  decay  time,  we  may  expect  a  suc¬ 
cession  of  (reflected)  waves  to  appear  in  the  west-central  basin,  with 
each  wave  less  intense  than  its  predecessor  and  gradually  merging  into 

the  western  boundary  current  ays tea.  The  corresponding  decay  time  with 

-1  -6  -1 
bottom  friction  is  t^  ■  k  ,  or  about  12  days  with  k  -  10  sec  In 

this  case  the  transient  waves  would  be  almost  completely  damped  in  a 

time  less  than  their  own  period,  and  would  effectively  disappear  from 

the  solution.  The  relative  prominence  (or  absence)  of  transient  Rossby 

waves  in  the  meridional  flow  in  the  western  part  ef  the  eceans  may 

therefore  be  indicative  of  the  relative  role  of  lateral  and  bottom 

friction,  although  a  nonlinear  viscosity  may  also  be  operative.  A 

similar  damping  ef  transient  motions  by  the  introduction  of  bottom 

friction  into  e  barocllnlc  oceanic  model  has  been  noted  by  Bryan  and 

Cox  (1968). 
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NONLINEAR  EFFECTS 

It  la  of  Interest  to  document  the  role  of  nonlinear  momentum  ad- 
vection  In  the  present  spherical  model,  and  to  this  end  a  comparative 
integration  was  made  with  the  omission  of  the  inertial  terms  in  (19) 
and  (20)  by  the  selection  S  -  0.  The  resulting  solutions  for  the  trans¬ 
port  at  30  days  are  shown  in  Fig.  8,  with  the  parameters  AX  -  A$  ■  2*, 
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A  -  10  cm  sec  ,  h  -  400  m,  and  with  the  zonal  stress  (26)  over  the 
60-deg  longitude,  48-deg  latitude  basin  as  before.  We  see  at  once  that 
the  transport  is  similsr  in  the  linear  and  nonlinear  solutions.  The 
maximum  eastward  flow  in  mid-latitudes  is  slightly  weaker  in  the  linear 
version,  but  elsewhere  the  zonal  transport  solutions  are  almost  indis¬ 
tinguishable.  The  maximum  meridional  boundary-current  transports  are 
likewise  very  nearly  the  same  in  both  solutions,  as  is  the  interior 
transport  pattern.  Only  the  meridional  countercurrents  in  the  west 
are  noticeably  weaker  in  the  linear  version,  and  this  effect  is  in 
accord  with  the  results  of  Veronls  (1966c)  for  the  corresponding  non- 
dlvergent  6-plane  model.  The  meridional  transient  motions  seaward  of 
these  countercurrents  are  likewise  weakened  by  the  deletion  of  the 
nonlinear  terms,  and  in  general  display  less  meridional  tilt  than  their 
nonlinear  counterparts.  This  effect  was  previously  noted  in  the  6- 
plane  solutions  discussed  by  Gates  (1968) ,  although  it  is  apparently 
less  marked  in  the  present  spherical  formulation. 

We  note  that  the  nonlinear  solutions  in  Fig.  8  are  for  the  value 
S  -  4/3,  which  corresponds  to  a  linear  decrease  of  speed  with  depth 
toward  zero  at  the  ocean  bottom.  The  solutions  with  this  value  do  not 
differ  significantly  from  those  with  S  -  1  (not  shown) ,  corresponding 
to  a  current  uniform  over  depth,  nor  from  those  with  S  -  4  (not  shown), 
corresponding  to  a  more  nearly  exponential  (and  possible  more  realistic) 
current  profile.  For  the  S  -  4  case  the  Interior  transient  motions  are 
slightly  less  intense  than  those  shown  in  Fig.  8,  although  the  overall 
similarity  of  the  solutions  is  testimony  to  the  relatively  minor  role 
played  by  the  nonlinearity  Introduced  by  S  of  the  erder  unity.  This 
is  in  accord  with  the  three-dimensional  solutions  of  Bryan  and  Cox 
(1967)  for  a  similarly  modest  degree  of  nonlinearity.  The  more  highly 
nonlinear  solutions  reported  by  Bryan  and  Cox  (1968)  would  appear  to 


5n  30  days 


30  da 


The  linear  solutions  for  zonal  transport  (a)  and  meridi- 
onsl  transport  (b)  In  Sv  (100  km)"1  st  30  days  with 
S  ■  0  (solid  llnss) ,  compared  with  the  corresponding  non- 
1  inear  solutions  with  S  -  4/3  (dashed  lines) .  Here  A  ■ 


10®  cm2  sec"1  and  h  -  400  m.  The  mean  currents  u/h  and 
v/h  In  cm  sec-1  are  given  by  multiplying  the  transport  lso- 
lines  by  the  factor  5/2. 
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correspond  to  values  of  S  of  the  order  of  10  to  100,  and  represent  a 
marked  confinement  of  the  (wind-driven)  transport  to  the  surface  layers. 
Integrations  for  this  range  of  S  have  not  been  performed  with  the  pre¬ 
sent  model. 

DEPTH  EFFECTS 

Since  the  ocean  depth  has  been  assumed  constant  in  the  present 
solutions,  its  value  has  been  selected  as  representing  in  some  way  the 
upper  oceanic  region,  wherein  the  wind-driven  barotroplc  mode  might  be 
expected  to  be  most  evident.  These  solutions  cannot  be  considered  ade¬ 
quate  simulations  of  the  oceanic  surface  layer,  however,  in  view  of 
their  neglect  of  the  substantial  variations  of  thermocllne  depth. 
Strictly  speaking,  we  have  here  studied  the  behavior  of  only  a  shallow 
homogeneous  sea,  and  it  may  therefore  be  of  Interest  to  examine  the  ef¬ 
fect  of  the  (uniform)  depth  selected  for  such  a  basin. 

The  comparative  solutions  for  zonal  and  meridional  transport  shown 

in  Fig.  9  display  the  effects  of  increasing  the  depth  h  from  400  m  to 
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1200  m,  with  the  values  S  ■  4/3  and  A  ■  2  *  10  cm  sec  common  to 
both  integrations.  The  solutions  for  zonal  transport  are  only  slightly 
altered  by  the  depth  Increase,  and  the  maximum  meridional  transports 
along  the  western  boundary  are  also  nearly  the  same  in  both  solutions. 
This  is  in  agreement  with  the  numerical  solutions  of  Bryan  and  Cox 
(1967) ,  which  show  the  mean  wind-driven  transport  to  be  relatively  in¬ 
sensitive  to  the  basin's  (uniform)  depth  provided  it  exceeds  a  charac¬ 
teristic  scale  depth.  There  are,  however,  noticeable  differences  in 
the  position  and  strength  of  the  transients,  which  are  most  clearly 
seen  in  the  west-central  basin  in  the  meridional  transport  solutions 
of  Fig.  9.  These  transients  are  Rossby  waves  reflected  from  the  west¬ 
ern  boundary,  and  are  propagating  westward  relative  to  the  zonal  cur¬ 
rent.  The  behavior  of  such  waves  is  essentially  described  by  the  linear 
lnvlscid  theory  for  the  basin's  free  oscillations,  as  shewn  for  example, 
by  the  numerical  studies  of  Gates  (1969).  The  primary  effect  of  a  depth 
Increase  on  these  waves  would  therefore  be  an  increased  (reflected) 
transient  wavelength  and  an  increased  westward  phase  speed,  both  Induced 
through  the  divergence  permitted  by  the  free-surface  variations.  The 


The  solutions  for  sonsl  transport  (s)  and  nsrldlonsl 
transport  (b)  In  Sv  (100  ka)"^  at  30  days  with  tha  depth 
h  ■  1200  a  (solid  lines)  end  h  ■  400  a  (dashed  lines) . 

The  corresponding  solutions  for  water  elevation  In  ca  ere 
also  shown  (c)t  with  the  contours  drawn  tven  20  ca  In 
both  cases.  Here  S  ■  4/3  and  A  ■  2  *  10®  ca*  sec”*.  The 
aean  currents  u/h  and  v/h  In  ca  sec-*  are  given  by  Multiply 
lng  the  solid  transport  Iso  lines  by  the  factor  S/6  and  the 
dashed  lsollnes  by  the  factor  S/2. 
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phase  differences  between  the  transients  in  ths  south-central  basin  In 
Fig-  9  at  30  days,  for  eusple,  are  consistent  with  the  phase  speed  in¬ 
crease  (from  1.3  ■  sec”*  to  1.5  ■  sec  *)  to  be  expected  from  the  depth 
change  froa  400  a  to  1200  a.  The  apparent  increaae  in  transient  wave¬ 
length  aay  also  be  accounted  for  by  reflection  froa  the  western  shore 
in  the  deepened  basin. 

Ve  aay  speculate  that  changes  in  the  transient's  behavior  similar 
to  those  in  Fig.  9  would  be  produced  if  the  basin's  depth  were  increas¬ 
ed  by  a  further  factor  of  three  (froa  1200  a  to  3600  a),  in  order  to 
approxiaate  the  average  depth  of  the  world's  oceans.  On  the  basis  of 
the  present  evidence,  however,  no  significant  change  would  be  expected 
in  the  transport  or  structure  of  the  laportant  western  boundary  cur¬ 
rents.  The  present  solutions  thus  provide  justification  for  the  con¬ 
tinued  use  of  a  relatively  shallow  basin  in  studies  of  a  berotropic 
ocean,  at  least  in  a  program  of  coaparative  integrations  for  tha  wind- 
driven  transport. 

In  terns  of  the  vertically  integrated  dynamical  systea,  (19)  to 
(21),  the  appearance  of  the  depth  h  in  the  nonlinear  inertial  tarao  ia 
a  relatively  uniaportant  dependence  (at  least  for  the  caae  of  unifora 
h  and  with  S  of  the  order  unity).  If  these  terns  ara  omitted,  we  note 
that  the  only  reaaining  explicit  appearance  of  tha  basin  depth  i a  aa  a 
coefficient  of  the  horixontal  pressure-force  terns.  Since  the  (verti¬ 
cally-integrated)  transports  u  and  v  are  esentially  unaltered  by  depth 
changes  in  the  present  solutions  (and  are  quasi -steady  in  tha  laportant 
western  boundary  region) ,  we  aay  conclude  that  the  magnitude  of  the 
free-surfaca  elopes  JC/JX  and  will  vary  inversely  with  changes 

of  the  (uniform)  depth  h.  This  la  confined  by  the  comparative  solu¬ 
tions  for  (  also  shown  in  Fig.  9  (lower),  where  contours  differing  by 
a  factor  of  three  nearly  coincide  in  the  400  a  mad  1200  a  caaoa,  while 
the  sero  contour  renal  ns  in  approximately  the  same  position.  We  aay 
further  note  that  the  water  elevation  along  tha  boundary  ia  also  re¬ 
duced  by  an  approxiaate  factor  of  three.  Ia  an  ocean  4  ka  deep  we 
night  therefore  expect  a  water  elevation  of  about  10  cm  across  tha  aala 
western  boundary  currents  with  the  present  wlad-etreaa  pattern. 
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VI.  STATISTICAL  PROPOtTIES  OF  TOE  SOLUTIONS 

MEAN  AW)  VARIABILITY  FIELDS 

While  the  present  Integrations  have  been  compared  on  a  synoptic 

basis  only  at  selected  tiati,  the  comparisons  are  generally  representa- 

e 

tive  after  an  initial  spin-up  period  of  about  25  days.  Since  in  moat 
caaea  the  integrations  have  been  carried  out  to  60  days,  the  transport 
and  elevation  solutions  have  been  averaged  over  the  period  30  to  60 
days,  and  the  corresponding  rms  (root-mean-square)  variabilities  deter¬ 
mined.  In  general  these  flelda  portray  the  differences  among  the  solu¬ 
tions'  structure  and  intensity  as  revealed  by  the  synoptic  comparisons 
discussed  above,  and  need  not  be  presented  in  detail.  It  la  of  inter¬ 
est,  however,  to  ehov  theee  atatlatics  for  at  leaat  a  reference  (spher¬ 
ical)  caae. 
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For  the  case  A  ■  10  cm  sec  ,  S  ■  4/3  and  h  ■  400  m,  the  distri¬ 
butions  of  the  30-  to  60-day  average  transport  and  water  elevation  and 
the  assecieted  rms  variabilities  are  shown  in  Fig.  10,  determined  from 
the  solutions  every  12  hours  during  thia  period.  Here  the  wind  stress 
le  given  by  (26),  and  the  60*-longltude  basin  extends  from  18  N  to  66  H 
with  e  2*  («  A*  •  AX)  separation  between  adjacent  transport  (or  eleva¬ 
tion)  grid  pointa.  The  general  similarity  between  the  average  fields 
and  the  instentaneous  solutions  for  this  caae  shown  previously  in  Figs. 
3,  4,  7,  and  8  is  testimony  to  the  near-steadiness  of  the  major  fea¬ 
tures  of  the  solutions  (and  to  ths  stability  of  ths  calculations). 
Seaward  from  the  western  boundary  sons  the  average  transport.!  are  the 
result  of  the  passage  of  the  successive  transient  Rossby  waves  which 
dominate  ths  solutions,  and  in  general  ere  in  an  approximate  geostrophic 
balance  with  the  average  surface  elevation. 

The  — —  rms  transport  (and  elevation)  variability  in  the  south¬ 
western  portion  of  the  basin  is  believed  to  be  caused  by  the  reflection 

*Thie  spin-up  is  defined  as  the  time  required  for  the  basin's 
tetal  kinetic  energy  to  first  decrease  after  the  major  western  bound¬ 
ary  currents  have  reached  maximum  strength  from  an  Initial  state  of 
rest  (see  Fig.  11). 


Fit.  10.  Tha  dlatrlbatioos  of  average  tonal  treaty ort  la  tv  (100  !■)“ 
(•),  meridional  trontport  la  Sv  (100  ka)’l  (b),  and  anrfaca 
a  lava  doo  la  ca  (r)  for  tba  period  30  to  40  day* ,  ebova  by 
tha  eoltd  laollnaa  for  a  rafaraaca  (epherlcel)  caaa.  Tba 
10-day  average  carraata  la  ca  eec“l  ara  given  by  aal tip lying 
tba  traaaport  laollnaa  by  tba  factor  S/2.  Tba  cecrmpeadlag 
raa  (root  neaa-egnara)  varlabllltlaa  daring  tba  peeled  30  to 
SO  daya  are  eboya  by  tba  daabad  llaaa  la  tba  mm  aalta. 

■ara  A  •  10*  car  etc  ,  S  •  4/3.  and  b  •  400  a. 
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from  the  western  (and  southern)  shore  of  the  transient  wave  energy  of 
the  aajor  aid-latitude  gyre.  The  role  of  the  orientation  of  the  west¬ 
ern  shore  in  such  a  reflection  process  has  recently  been  examined  by 
Cates  (1969),  who  shows  that  the  transient  variability  nay  be  localized 
in  different  regions  of  the  (western)  basin  by  an  oblique  western  bound¬ 
ary.  In  general  the  ms  variability  of  the  aeridlonal  flow  is  approxi- 
aately  twice  that  of  the  zonal  notion,  and  is  coaparable  with  the  aver¬ 
age  aeridlonal  notion  itself  over  auch  of  the  basin  outside  the  western 
countercurrents.  The  me  variability  of  the  surface  elevation  given  in 
Mg.  10  is  coapatible  with  the  variability  of  the  (geostrophic)  aerid¬ 
lonal  transport  end  the  characteristic  (reflected)  transient  wavelength 
of  about  1600  ka  noted  earlier  for  this  case. 

In  the  dlecusalon  of  frictional  effects  (see  Pig.  7),  it  was  noted 
that  the  addition  of  bottoa  friction  to  the  lateral  viscous  nodal  sig¬ 
nificantly  daaped  the  aaplitude  of  the  transient  Rossby  waves,  end  re¬ 
duced  the  aegnltude  of  the  western  boundary  transports.  The  rae  vari¬ 
abilities  in  this  case  are  soaa  10  to  20  percent  of  the  values  shown  in 
Fig.  10.  In  the  cooperative  6-plane  integrations  (Figs.  3  and  4)  the 
rae  variabilities  ere  systeaatically  larger  than  thoee  shown  here  in 
the  northern  part  of  the  basin,  an  effect  caused  by  the  6-plene  geo- 
aatry .  The  variabilities  in  the  reaalnlng  cooperative  integrations 
discussed  earlier  (for  the  sane  basin  dlaenslons)  ere  all  coaparable 
to  thoee  shown  for  the  reference  case  in  Fig.  10. 

BWBCT  PAKTITIOWmC 

The  total  kinetic  energy  of  the  horizontal  flow  is  e  useful  mee- 
svre  of  the  behavior  of  a  wind-driven  ocean,  particularly  to  depict 
the  approach  of  the  dominant  boundary  currents  toward  e  steady  state 
(see,  for  exaaple,  Bryan  (1963),  or  Veronle  (1966b,  1966c).  In  Fig. 

11  the  total  kinetic  ead  potential  energies  ere  shown  as  e  function 
of  tlaa  for  the  reference  (spherical)  case;  they  were  computed  from 
the  expressions 


Kinetic 


Tim*  in  tfop- 


11.  The  variation  of  tha  total  kioatic  and  potential  energy  foi 
the  reference  caae  of  A  ■  10®  cm*  aec~l,  S  •  4/3,  h  •  400  ■ 
and  AX  •  &4  •  2*.  The  inaat  shove  the  typical  fine-scale 
behavior  in  which  the  kinetic  eoergy  ia  plotted  every  2  hr. 
The  daahed  curve  la  the  kioatic  energy  for  the  caae  of  A  • 
2  ■  10®  c*2  eec“l. 
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2  A  \ 

A  p  AX  m  %  a  2 

total  kinetic  energy  -  — — £^(sln  ♦j+i  "  sln  *j-i)(uij  +  vij)  <32> 


2 

e  gp  AX  _ 

total  potential  energy  -  - j - 2^  (»in  *j+1  ’  sln  ^j-i)  (33> 


where  the  sunmatlons  extend  over  all  transport  (u,v)  and  elevation  (O 
points  respectively.  Here  a  spln-up  or  Initial  energy  growth  period 
of  some  24  to  26  days  occurs,  after  which  the  energy  shows  damped  oscil¬ 
lations  with  free  period  of  the  baaln  (approximately  13  days).  Both 
the  kinetic  and  potential  energy  characteristically  fall  to  achieve 
their  maximum  until  the  second  (free)  period.  This  behavior  Is  In  con¬ 
trast  to  that  of  the  6-plane  solutions,  which  typically  reach  maximum 
energy  during  the  first  period  (see,  for  example,  Cates  [1968]).  This 
deferred  spln-up  also  occurs  In  the  caae  of  increased  lateral  viscosity. 
Superposed  on  theae  curves  Is  a  fine-scale  oscillation  illustrated  by 
the  lnaet  In  Pig.  11,  where  a  period  of  approximately  13  hr  may  be  iden¬ 
tified  with  the  passage  of  (long)  surface  gravity  waves  across  the 
basis.  The  rms  currents  associated  with  these  waves  are  evidently  of 
the  order  of  0.02  cm  sec  \  and  are  much  lees  than  the  typical  currents 
of  about  1  cm  sec  1  aseociated  with  the  transient  Rossby  waves. 

We  may  also  note  that  the  potential  energy  is  here  approximately 
half  the  kinetic  energy,  by  virtue  of  the  relatively  shallow  basin 
depth  selected.  Recalling  that  the  trsnsporte  u  snd  v  are  nearly  In¬ 
dependent  of  variations  of  the  basin's  (uniform)  depth  while  the  surface 
elevation  varies  Inversely  with  h  (see  Fig.  9),  the  total  kinetic  energy 

(32)  will  vary  as  h”*  and  the  total  potential  energy  (33)  will  vary 
-2 

as  h  ,  other  factors  being  the  came.  For  a  basin  of  depth  h  •  4  km  we 
could  therefore  expect  the  potential  and  kinetic  energies  to  be  in  the 
ratio  of  about  1/20,  rather  than  the  ratio  of  1/2  In  the  present  model. 

Although  the  bulk  of  the  kinetic  energy  le  evidently  associated 
with  the  western  boundary  currents,  it  is  of  Interest  to  make  a  more 
quantitative  assessment  and  to  examine  the  energy  of  the  transient 
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waves  in  particular.  In  TabLa  4  the  partitioning  of  tha  kinetic  energy 
of  both  zonal  and  meridional  aotlona  among  a  zonal  nean,  standing  zonal 
eddies,  and  transient  zonal  eddlae  is  shewn  for  the  reference  solution 
(see  Pigs.  10  and  11).  Hera  the  zonal  average  is  denoted  by  the  opar- 
etor  [(  )],  the  departure  fron  the  zonal  average  by  (  )*,  the  tine  av¬ 
erage  (determined  fron  the  solutions  avery  12  hr  during  the  period 
30  to  60  days)  by  (  ) ,  and  the  departure  fron  the  tine  average  by  (  )'. 

In  the  basin  as  a  whole  approximately  71  percent  of  the  total  kinetic 
energy  is  seen  to  reside  in  the  nerldlonal  notion  of  the  standing  waves, 
principally  the  western  boundery  current.  The  nean  zonal  notion  and 
the  zonal  notion  of  the  standing  waves  (which  complete  the  major  gyral 
circulation  of  the  basin)  account  for  about  28  percent  of  tha  total  ki¬ 
netic  energy.  The  tranelart  waves  (here  the  zonally  moving  kossby 
waves)  represent  only  about  1  percent  of  the  total  kinetic  energy, 
which  is  Itself  apportioned  between  meridional  and  tonal  notion  in  ap¬ 
proximately  the  ratio  3:1.  In  spite  of  this  dominance  of  the  total 
energy  by  the  steady  nodes  (the  potential  energy  is  similarly  appor¬ 
tioned),  the  relatively  weak  tranaients  nevertheless  control  the  meri¬ 
dional  transport  over  most  of  the  area  of  the  basin,  as  seen,  for  ex¬ 
ample,  in  Pig.  3.  The  dsta  of  Table  4  also  show  that  the  transient 
waves'  contribution  to  the  kinetic  energy  is  greatest  in  the  southern 
portion  of  tha  basin,  a  feature  which  was  also  noted  previously. 

AHCULAK  HOHPfTUH  FLUX 

further  Insight  into  the  role  of  transient  planetary  waves  in  a 
wind-driven  circulation  is  afforded  by  the  processes  which  are  respon¬ 
sible  for  the  flux  of  angular  momsntun  in  the  basin.  Table  S  shows 
the  poleward  flux  of  tonal  angular  momsntun  during  the  period  30  to  60 
dsye  for  the  reference  case  of  figs.  10  end  11  and  Table  4,  partitioned 
among  the  fluxes  due  to  the  tonal  nean  motions,  the  standing  tonal  waves, 
and  the  transient  tonal  waves.  The  standing  eddies  are  seen  te  sc count 
for  about  95  percent  of  the  total  mo  man  turn  transport,  and  have  s  mnimrn 
flux  convergence  near  the  latitudes  of  ths  maxims  of  the  western  bound¬ 
ary  currents  (see  fig.  10);  the  bulk  of  this  stationary -eddy  flux  con¬ 
vergence  nay  therefore  be  assoc is ted  with  the  western  boundary  current 
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system.  The  flux  of  the  mean  tonally  averaged  motions  [u][v]  is  rela¬ 
tively  small  because  of  the  smallness  of  [v] ,  which  is  in  turn  a  reflec¬ 
tion  of  the  near  steadiness  of  the  30-  to  60-day  mean  flow.  The  tran¬ 
sient  (Rossby)  waves  produce  a  northward  (positive)  momentum  flux  over 
the  latitudes  occupied  by  the  southern  or  principal  gyre,  with  a  negative 
flux  in  the  northern  gyre.  This  flux  distribution  is  a  reflection  of 
the  systematic  meridional  tilt  of  the  transient  waves  during  this  per¬ 
iod,  which  is  characteristic  of  all  the  present  solutions  for  basins 
with  a  meridional  western  boundary  (see,  for  example.  Fig.  3).  The 
convergence  of  this  transient-eddy  flux  is  a  maximum  near  the  gyral 
boundary,  where  the  mean  tonal  stress  [t^]  is  a  maximum  and  where  the 
aman  tonal  transport  [u]  is  strongest. 

On  the  basis  of  these  (end  elmllar)  calculations,  the  transient 
waves  may  be  said  to  contribute  in  a  small  but  systesuitlc  way  to  the 
maintenance  of  ths  mean  tonal  circulation  In  the  basin,  in  the  same 
sense  that  the  atmospheric  Rossby  waves  contribute  to  the  maintenance 
of  the  mean  tonal  winds.  In  the  atmospheric  case,  however,  the  tran¬ 
sient  eddy  flux  Is  somewhat  larger  than  the  flux  of  the  tonal  mean  and 
standing  eddies,  and  the  flux  divergence  is  s  good  measure  of  the  mean 
tonal  stress  of  the  atmosphere  on  the  earth  (sse,  for  example,  Starr 
(1968)  or  Holopainen  [1967]).  In  the  ocean  (or  more  precisely,  in  the 
present  ocean  model),  on  the  other  hand,  the  standing  eddy  flux  is 
much  larger  than  the  flux  of  the  other  modss.  Ths  present  oceanic 
case  la  further  distinguished  from  the  atmospheric  regime  by  the  fact 
that  the  total  oceanic  flux  convergence  by  all  nodes  (see  Table  5)  is 
apparently  not  an  important  process  in  the  maintenance  of  the  tonal 
momentum  balance.  Here  the  applied  zonal  wind  stress  t*  (or  [t^]  is 
approxlsmtely  balanced  by  a  net  pressure  torque  exsrted  against  the 
eastern  and  western  walls  of  the  basin  by  ths  wind-raised  water,  with 
the  effects  of  latsrsl  friction,  mean  meridional  motion,  end  meridional 
flux  convergence  all  sawller  by  several  orders  of  magnitude. 

Although  s  couplets  study  of  the  kinetic-energy  belance  has  not 
been  made,  the  eddy  momentum  flux  discussed  above  is  associated  with 
the  transformation  between  the  kinetic  energy  of  the  treneient  eddies 
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and  Chat  of  Che  mean  zonal  aotlon.  If  wa  assume  that  this  energy  trans¬ 
formation  la  proportional  Co  [u' v' ]  3[u]/3y,  by  applying  Che  enalyala 
of  UcbsCer  (1961),  for  exanple,  Co  Che  present  vertically  Integrated 
aodel,  ve  note  froa  Che  data  of  Table  5  Chat  this  product  is  positive 
across  virtually  the  entire  latitudinal  width  of  the  baaln.  At  leeat 
on  the  zonal  average,  therefore,  the  transient  Rossby  waves  may  be  re¬ 
garded  as  sys teas t leal ly  supplying  energy  to  the  aean  zonal  aotlon. 

As  was  the  case  in  the  aoaentua  balance,  however,  this  eddy  contribu¬ 
tion  Is  only  a  saall  part  of  the  kinetic  energy  budget  of  the  aeen 
flow;  here  the  doalnant  effects  are  the  direct  generetlon  of  the  kinet¬ 
ic  energy  of  the  aean  (zonal)  notion  by  the  zonal  wind  stresa  and  lta 
dlsalpatlon  by  the  lateral  viscosity,  which  occurs  principally  in  the 
standing  eddy  which  is  the  western  boundary  current  Itself.  The  pre¬ 
sence  of  nerldlonal  boundaries  thus  aakes  the  character  of  the  wind- 
driven  oceanic  circulation  quite  different  froa  the  corresponding  cir¬ 
culation  in  a  laterally  unbounded  (global)  ocean,  or  froa  the  analogous 
barotroplc  etaoapherlc  circulation  reglae.  The  large  zonal  boundary 
torque  and  the  proalnent  standing  (western)  boundary  currents,  with 
their  dominance  of  both  the  kinetic  energy  end  aoaontaa  budgets,  would 
presumably  be  absent  in  a  zonal  (circumpolar)  ocean  channel.  The  role 
of  the  transient  planetary  waves  In  the  maintenance  of  the  Antarctic 
circumpolar  circulation,  for  exanple,  might  therefore  be  expected  to 
be  more  prominent  than  it  le  In  the  tonally  bouoded  major  ocoaa  baa  ins . 
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VII.  CONCLUSIONS 

Th«  comparative  turner leal  integration*  preaented  here  aerve  aa  a 
preliminary  documentation  of  the  behavior  of  a  wind-driven  homogeneous 
ocean  In  a  basin  of  uniform  depth  on  the  sphere.  The  examination  of 
the  model's  sensitivity  to  changes  of  the  finite-difference  grid  reso¬ 
lution,  to  the  use  of  6-plane  geometry,  and  to  the  basin  size  serve  to 
Illustrate  the  character  of  the  purely  numerical  or  geometrical  effects 
Inherent  In  such  a  model.  From  theae  teats  we  may  conclude  that  at 
least  a  2-deg-latltude  grid-mash  resolution  la  required  to  portray  satis¬ 
factorily  the  major  current  features;  that  the  B-plane  approximation 
la  adequate  for  the  middle  and  lower  latitudes;  and  that  the  western 
boundary  current  transport  la  approximately  proportional  to  the  tonal 
width  of  the  basin.  The  model's  response  to  changes  of  Its  physical 
parameters  representing  the  effects  of  friction,  nonlinearity,  and 
water  depth  has  also  been  examined.  From  these  tests  we  may  conclude 
that  an  Increase  In  the  magnitude  of  the  lateral  eddy  viscosity  gener¬ 
ally  reduces  the  strength  of  the  western  boundary  currents;  that  the 
introduction  of  bottom  friction  reduces  the  Intensity  of  both  the  west¬ 
ern  boundary  currents  and  the  transients;  that  the  vertically  Integrated 
transport  Is  only  slightly  affected  by  the  inertial  nonlinearity  repre¬ 
sented  by  simple  assumed  current  profiles;  and  that  sn  Increase  in  the 
basin's  (uniform)  depth  has  little  effect  upon  the  total  transport,  but 
produces  an  approximately  proportional  decrease  of  the  free  surface 
water  elevation  and  slope. 

These  conclusions  apply  to  the  quasi-steady  circulation  produced 
by  a  large-scale  tonal  wind  stress  In  s  basin  with  purely  meridional 
and  tonal  boundaries.  The  transient  plsrotary  circulations  excited  by 
this  stress  are  the  low-order  Rossby  wave  nodes  and  the  transients  prom¬ 
inent  In  the  west-central  basin  after  about  20  days  may  be  Identified 
aa  the  reflections  of  these  waves  from  the  western  shore  (Cates,  196S). 
Theae  waves  are  responsible  for  a  transport  variation  in  which  the  rms 
variability  of  the  meridional  current  Is  systematically  larger  than 
the  mean  meridional  flow  over  a  considerable  portion  of  the  western 
basin,  and  la  usually  larger  than  the  rms  variability  of  the  tonal 


-51- 


curreat  es  well.  Changes  of  the  bee  In  site  or  depth,  or  aodest  Increases 
In  the  saount  of  leterel  friction  or  inertial  nonlinearity,  have  rela¬ 
tively  little  effect  upon  the  transient  waves  in  the  central  basin. 

This  is  in  accordance  with  the  identification  of  these  notions  with 
free  Rossby  waves,  whose  structure  and  behavior  are  described  with  good 
approximation  by  the  linear  inviscld  equations.  The  amplitude  of  the 
transients,  however,  is  markedly  reduced  by  the  Introduction  of  bottom 
(Ekaan)  friction. 

In  ell  of  the  present  integrations  the  basin  has  been  assumed  to 
be  of  uniform  depth,  end  the  wind  stress  has  been  assumed  steady  (aa 
well  as  tonal).  The  Introduction  of  time -dependent  stress  nay  be  ex¬ 
pected  to  enhance  the  role  of  the  transient  wavee  in  the  maintenance 
of  the  naan  circulation  through  resonance  and  nonlinear  effects  (Ped- 
losky,  1965,  1967;  Veronia,  1966a),  snd  tha  presence  of  bottom  topo¬ 
graphy  may  likewise  exert  an  important  influence  on  both  the  mean  and 
transient  notion  in  a  homogeneous  ocean  (Holland,  1967).  The  resulta 
of  numerical  Investigations  of  these  effects  will  be  reported  sepa¬ 
rately. 
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